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I. INTRODUCTION 



Recently, we have calculated the complete set of massive one-loop master integrals 
needed in the calculation of the next-to-next-to-leading order (NNLO) parton 
model corrections to the hadropro duct ion of heavy flavors 0. We used Feynman 
parametrization to evaluate the one-loop master integrals in n = 4 — 2e dimensions. 
We obtained the coefficients of the Laurent series expansion of the relevant scalar 
integrals in terms of the parameter e up to 0(e 2 ) as needed for the NNLO calculation. 
We found that the real parts of some of the e 2 coefficients contain a new class of 
functions which can be written in terms of one-dimensional integral representations 
involving products of log and dilog functions. These so-called single and triple index 
L functions cannot be expressed in terms of classical polylogaritms but can be seen 
to belong to a generalization of the classical polylogarithms which are called multiple 
polylogarithms. 

Functions analogous to the triple index functions L ai(T2aa also arise in the ap- 
proach of [3j when one analytically continues their 0(e 2 ) integral representation for 
a general vertex function. Methods differing from ours have been used for the deriva- 
tion of master iV-point integrals such as the differential equations method |I] or the 
nested sum method |5... Depending on the number of scales involved, the results in- 
clude multiple polylogarithms 6J and/or harmonic [7j or two-dimensional harmonic 
jH] polylogarithms. The latter function all are subsets of multiple polylogarithms. 
Presenting our results in terms of multiple polylogarithms will facilitate a compar- 
ison with the results of possible rederivations of the scalar one-loop integrals using 
other methods. It is very likely that future results of multiloop calculations will be 
presented in terms of multiple polylogarithms or their subclasses. Alongside with 
this the necessary tools will be developed to deal with multiple polylogarithms, be 
it analytically or numerically. In fact, recently a computer code has been written 
for the numerical evaluation of the multiple polylogarithms jHj- It is therefore timely 
that we express the results of PQ also in terms of multiple polylogarithms. 

It is a purpose of this paper to show that the single and triple index L functions 
introduced in [1 can all be related to multiple polylogarithms. This is done in explicit 
form. We are thus able to present our results for the scalar massive one-loop master 
integrals in terms of multiple polylogarithms and classical polylogarithms j^j. In 
Sec. II we recapitulate material on the definition of the single and triple index L 
functions as they arise in the approach of pQ. Simple symmetry relations allow one 
to restrict the discussion to the triple index L functions and L +++ , and to the 

single index L function L + . In Sec. II we also recapitulate the definition of multiple 
polylogarithms. In the subsequent sections we will write down the formulas needed 
to transform the L functions to multiple polylogarithms for general arguments. The 
general formulas are not always applicable when the arguments take special values 
as they do in the massive one-loop calculation. For these special values one must 
carefully discuss the limiting behavior of the general formulas. In Sec. Ill A we 
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derive the general formula which relates the L h+ functions to the set of multiple 

polylogarithms. Section III B considers special cases of the general relation. Similarly, 
Sec. IV A gives general relations which allow one to express the L +++ functions 
in terms of multiple polylogarithms. In Sec. IV B we discuss special cases for the 
arguments of the L +++ functions. Sections V A and V B repeat the discussion for 
the single index L + functions. Finally, Sec. VI presents our conclusions. 

As remarked on before, the L functions appear only in the real parts of some of 
the 0(e 2 ) coefficient functions of the masive one-loop integrals. In the notation of PQ 

(2) (2) (2) 

these are the three-point coefficient functions ReC\ , ReC^ and ReC^ , and the 

(2) (2) (2) 

four-point coefficient functions Re-D^ , Re_D 2 and Re-D^ . For the sake of brevity 
we have decided to present multiple polylogarithm results in this paper only for the 
four-point coefficient function ReD[ . This result is listed in the Appendix. The 
corresponding results for the other five coefficient functions are readily available in 
electronic form 



II. BASIC FEATURES 

In order to make the paper self-contained, we write down a number of basic 
definitions for the L functions and the multiple polylogarithms in this section, as 
well as some symmetry properties and domains of definitions for the single and triple 
index L functions. These will be of help when presenting the subsequent material. 

The definition for the L functions is as follows pQ: 

f 1 ln(«i + a x y) ln(a 2 + o 2 y) ln(a 3 + a 3 y) 
L ai(T2a3 [ai, a 2 , a 3 , a A ) = / dy ■ (1) 

Jo ct4 + y 

and 

T i N f 1 \n(ax + a 1 y)U 2 (a 2 + a 3 y) 

L <Tl {a 1 ,a2,a 3 ,a 4: ) = dy . (2) 

Jo «4 + y 

Here the (i = 1, 2, 3) take the values ±1 and the <x,-'s are either integers {1, 0, —1} 
or else kinematical variables. We want to emphasize that the numerical evaluation of 
the L functions is straightforward. 

The L functions possess simple symmetry properties as follows. One notices that 
a change of the integration variable y — > 1 — y results in the identity 

L ai (a 1} a 2) a 3) a i ) = -L_ ai (aci + cr 1 ,a 2 + a 3 , -a 3 , -a 4 - 1) (3) 

which implies that L_ can always be related to L + , and vice versa. We have thus 
written our results for the three-point and four-point functions in PQ only in terms of 
the L + functions. 

Turning to the triple index L function one notices that L UlU2aa (ai, a 2 , a 3 , a^) is 
symmetric under permutations of any two pairs of indices and arguments {<7j, a{\ and 
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{(Xj, otj} for (i 7^ j). The same change of variables as above y — > 1 — y results in 

£<TKT2CT 3 (a!i, "2, «3, "4) = —L- ai - a2 - a3 (ai + cti, ct 2 + cr 2 , ct 3 + cr 3 , — a 4 - 1). (4) 

Therefore, from the eight functions L , L h , L H , L_++, £++-> 

and only two are independent. We have chosen to write our results in terms of 

and L +++ . 

The domains of definition of the functions and L + that follow from 

the requirement that these functions take real values can be read off from the defining 
relations Eqs. (JTJ) and (J2J) considering the arguments of the log and dilog functions 
in the integrands, as well as from ensuring that the denominator of Eqs. (pp) and (j2j) 
does not change sign on the integration path. One has 



L+++(dll, 0>2, "3, "4; 

L-++(oix, "2, a 3 , 014) 

L+(a>i, a>2, "3, "4y 



ax > 0, «2 > 0, 0:3 > 0, a 4 < — 1 or 0:4 > 0; 
ax > 1, a2 > 0, 0:3 > 0, 0:4 < — 1 or 0:4 > 0; (5) 
cti > 0, «2 < 1, 0:2 + a 3 < 1, 03 7^ 0, «4 < —1 or 04 > 0. 

Looking at the definition of the triple index L function in (Q) one concludes that the 
boundary points a% — and/or a 2 = and/or a 3 = can be included in the domain 
of the definition for L +++ . The same holds true for a.\ = 1 and/or a 2 = and/or 
a 3 = for L_ ++ . Also, from the definition of the single index function L + in (j2J) one 
concludes that the boundary point ai = can be added to its domain of definition. 

The points 0:4 = { — 1,0} can also be included in the domain if the values taken 
by the other parameters guarantee the convergence of the integral. We mention 
that for all of our purposes the conditions with the boundary points included, 
are satisfied, e.g., our results for the integrals are real. Nevertheless, it is of course 
always possible to analytically continue the parameters to the complex plane. 

There are some further relations for the L functions which result from applying 
integration-by-parts identities. They are not listed here but can be found in Ap- 
pendix C of p. They have been used to reduce the set of L functions occurring in 
the master integrals to a subset of L functions having real values in physical phase 
space jl]. 

Multiple polylogarithms are defined as a limit of Z sums 0, e.g., 



Lim k ,...,m 1 (Xk } X\) — lhn ^2 mi m 2 m k • (6) 



m>n 2 ...>n fe >0 ,h "" ) ■■■' Ll 



The number w = mi + ... + is called the weight and k is called the depth of the 
multiple polylogarithm. The power series (jUJ) is convergent for \xi\ < 1, and can be 
analytically continued via the iterated integral representation: 



xiX2---%k / , \ mi— 1 j. 



-O 







t J X2X 3 ...X k -t 

(7) 



'dt \ m2_1 dt (dt \ m ^ 1 dt 

— o o ... o — o , 

t J x 3 ...x k -t \t I l-t 
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where the following notation is used for the iterated integrals: 

dt dt r dt n f dt n -i f dti 



o ... o 



(| a n — t <x\ — t 7 a n — t n J a n -i — t n _i j a% — ti 



x ... x (8) 



III. TRANSFORMATION OF L _ ++ TO MULTIPLE POLYLOGARITHMS 

In this section we will show that all our L_ ++ functions can be expressed in terms 
of multiple polylogarithms. 



A. General case for the L_ ++ function 
We begin with the function Eq. (JTJ), 

T , \ ) a ln ( tti ~ y) ln ( a2 + y) ln (" 3 + y) m\ 

L- ++ (ai,a2,a 3 ,a 4: ) = / dy ■ . (9) 

o 

After changing the integration variable y = ait one gets 

l/a\ l/«i 

^ln(«i — a±t) ln(a 2 + oc\t) ln(a 3 + ait) r ^ln«i ln(a 2 + ait) ln(«3 + ait) 

24 +t ~~ J 2i +t 

Y Ml - Wnai + ln(g + t)][lnai + ln(g + t )] = 

J 2i +t 

n "i 



InaJ V^ + ^^ + ^ lnV / Vf " * } 

V r ai ln(l-mn(2i + t) ^ ln(l - t) ln(23 + t) 

+ lna 1 / dt ' ^ V ina, / dt (10) 

ai o ai 



2i +t 

a i 

With the help of (J7|) the integral in the second term of the last equation of (|lUj) can 
be written as 

V v^= T^— / T ^-=Li 1 i(-^ --). ai) 

2, + t 7 l"*a 1,U a x ' a 4 j 1 ] 
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The third and fourth terms of the last equation in (fTU|) contain integrals of the form 

y<i-«)Nft+o, (12) 

To express such integrals in terms of multiple polylogarithms one proceeds as follows: 



fa Pi J Jfa + hJ p 2 + ti J p 2 + h J f3x + t 2 

U U U ti 

ln(/3i + t m ) / dtx—— / dti -—— = (13) 

-t m \ V ln(l-t 1 )K3i + *i) 



\n(Pi + t m )Lii,i - J dh 



In the first line of (|13|) we have changed the order of integration in the two-dimensional 
integral. We shall frequently use this trick further on. From Eq. one immediately 
concludes that 

dt — = L, w f -ft, -, — 1 + In (A + t m ) L. lfl f-A, ) ■ 

u 

(14) 

Let us now turn to the more involved integral [first term of Eq. (fTUj) ]: 

1 — 1/Q2 

d lnjaz + y) \n(a 3 + y) y ^-_ a2 t _ j ^ [In a 2 + ln(l - i)] ln(a 3 - a 2 i) _ 



, , ln(a 3 ~ a,t) ln(l - t)[lna 2 + ln(^ - t)] 
-lna 2 y dt ^— J dt ir — a = (15) 

+ ,„ a2 / ^Mas + rf _ , n a2 T d M^A _ 7 V (1 " Mg - 



a2 02 

The integral in the first term can be expressed as 

1 - 1/(>3 

ln(q 3 + y) j/^-a 3 i _ /" ^J lna 3 + ln(l - t)) _ 



Inaalnf^+L^f^-J-). (16) 
V a 4 / Va; 3 a 3 / 
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The integral in the second term can be written as 

— l/aa 



The third term from the last line of Eq. (fTK|) has a form which is an analogue of the 
integral (JT2j) and can be calculated in a similar way, 



f jn(l-t) Wjh-j) ( fa t m \ I t m \ 

Combining the Eqs. (fIH)l. (fTTj) and (fTH|) we arrive at the result for Eq. (fTHjl . 

1 ln(a 2 + y) ln(a 3 + y) /a 4 a 3 1\ . /a 4 1 
ay ■ = — , — , + lna 2 ^i,i — , 



a 4 + y ' ' \a 2 a 4 a 3 / ' \a: 3 a 4 

u 

+ ln(l + a 3 )^i,i f — , -—) + lna 2 lna 3 In (^±±1) . (19) 

\«2 «4/ V « 4 / 

Because the initial integrand is symmetric under the exchange of the parameters a 2 
and a 3 , the rhs of (frT^) can be rewritten in a symmetric form if desired. 

We are now left with the fifth term in ()1()J) . The fifth term is an integral of the 
type 

tm dt ln(l - t) ln( 7 i + t) ln( 72 + Q , . 

73 + t ' 1 J 

u 

In order to express such integrals in terms of multiple polylogarithms one can perform 
the following chain of transformations resulting in a multiple polylogarithm of weight 
four: 

T . ( 72 7i -t m \ y dt 4 J dt 3 f dt 2 J dti . . 
- ^1,1,1,1 -7s, — , — , = / — / — / — / i = 21 

V 73 72 71 / J 7l + U J 72 + *3 J 73 + *2 J 1 - *1 

dh f dt 3 J ln(l - 1 2 ) = _ *? dt 4 J d ln(l - 1 2 ) J dt 3 
7i + t 4 / 72 + £3 / 2 73 + ^2 / 7i + *4 / 2 73 + *2 / 72 + t 3 

*2 

f[ 7l + ^4 jj 73 + ^2 J 7l + ^4 J 73 + ^2 

+ t L t2 Ml-t 2 )ln( 72 + t 2 ) = 
J 13 + h J 7i + U 

12 

ln( 7 i + t 2 ) ln( 72 + t 2 ) ln(l - t 2 ) 



-I'(t m ) + I"{t m ) - J dt 2 



73 + ^2 
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where we have introduced the notation 







J"(t m ) = ln( 7l + t m ) J dt 2 



ln(l - t 2 ) ln(7 2 + h) 



n 73 + ^2 



The third term on the last line of (}2T|) is exactly the integral of the required type 
Eq. (jUJ). 

The integral in I"(t m ) has the form of (|14jl. For the integral I'(t m ) we write 

r (U = hu^^A h t2 ^pl = 'j dt ^2±lA Ul , U 3 ,=h) . (23 ) 

J 7l + *4 { 73 + ^2 J li+h V 73 / 

On the other, hand one has 

^1,1,1,1 -73, — , — , = / — / -73, = 24 

V 73 7i 72 / J 72 + h J 7i + h \ 73 / 

dti T . ( —ti\ 7 dt 2 



^ 7i + *i V 73 / ^ 72 + h 



ln(7 2 + t m ) / dt \ Li 1 1 (-73, — ) - / dt j—^ + p- Lii A (-73,— 
^7i+^i \ 7s / J 71 + *i V 73 , 



- ln(7 2 + t m )Li 1AA ( -73, — , — ^ ) - /'(*, 
V 7s 7i/ 



One then concludes that 



= -^1,1,1,1 ( -73, — , — , — ^ J - hi(7 2 + t m )Lii X i I -73, — , — ^ ) . (25) 
V 73 7l 72 / V 73 7i / 

Finally, substituting I'(t m ) and I"(t m ) into Eq. (|21|) we write down the result for the 
integral of the required type Eq. (J2UJ) . 

-m ^ ln(l-t) Hrfi+t) ln( 72 + t) . , f ~M 
a* — : = hi(7i + t m ) ln(7 2 + t m )Ln t i -73, + 

73 + * V 73 / 

ln(7 2 + t m )L«! 1 1 ( -73, — , ) + ln(7i + £ TO )Lii,i,i ( -73, — , — ^ J (26) 
V 7s 7i / V 7s 72 / 



. T ■ ( 72 7i -*m \ . T . ( 71 72 -t 

+-^1,1,1,1 -73, — , — , + 1,1,1 -73, — , — 



73 72 7i / \ 73 7i 72 
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We are now in the position to collect all required contributions to express the 
function in terms of multiple polylogarithms. Taking into account Eqs. (fTTjl. 
(JUJ), (Jinj), and (J23) we obtain 

/ a 4 a 2 «3 1 

L_ ++ («l, tt 2 , «3, «4j = ^1,1,1,1 , , , 

V ai a 4 a 2 «3 

/ OJ4 OJ3 «2 1 \ r / «4 «3 1 

> — > — > + lnaiLzi n — , — , 

V a\ «4 «3 d 2 J V«2 «4 «3 

+ ln(l + a 2 )Lii,i,i (-— , — , — + ln(l + a 3 )Lii,i,i (-— , — , — - 
V o.\ «4 0:3/ V ai a 4 «2 

+ lnai lna 2 Lii 1 (—, — — ) + lnai ln(l + a 3 )Lii 1 (—, — — ) (27) 
Va 3 a 4 / \a 2 «4> 

+ ln(l + a 2 ) hi(l + 0:3)^1 1 f 1 + In oc\ In «2 In 03 In ( 4 

' V «i 04/ V «4 

Some remarks are in order at this place. The final formula (|2*7|) contains multiple 
polylogarithms up to weight four. All multiple polylogarithms up to weight three 
can be expressed in terms of logarithms and clasical polylogarithms Li 2 and Li 3 . 
This fact is used by us when we reexpress our results for the massive scalar integrals 
in terms of multiple polylogarithms, i.e., our final results will contain only multiple 
polylogarithms of weight four. For the variables ati the conditions © are assumed. 
But in the results for the massive scalar integrals there are also cases when a% = 1 
and/or a 2 = and/or 0J3 = and/or a 4 = {—1,0}. In such cases the general formula 
(|27jl is no longer valid and these cases must be studied separately. 



B. Special cases for the L _ ++ function 

In the Laurent series expansion of the massive scalar one-loop integrals one en- 
counters special values of the arguments a, for which the general formula Eq. (|27j) 
no longer applies. This is quite obvious from the list of special cases discussed in the 
following. 



1. a\ = 1, CK4 = 

In such case one can make use of Eq. (|26|). One should find the limit of the 
expression on the right-hand side for t m — 1,73 — > 0. One obtains 



1 

/ 



dt 



ln(l - t) ln( 7l + t) ln( 72 + t) _ 
t 

dt 2 f dti r dt$ f dt 2 f dt\ 



Mm { ln( 7l + 1) 1„( 72 + 1) / / + 1„( 72 + 1) / I I , 

73^0 I J -73 - t 2 J 1 - ti J q -71 - t 3 J -73 - t 2 J Q 1 - h 
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^3 ^2 1 ^4 ^3 ^2 

o?t 3 /" ^2 /" dti f dt4 f dt 3 f dt2 f dt\ 



+ r dt 3 r dt 2 r^ti_ + r d U j dt 3 r 

j -72 -hJ Q -73 - h J Q 1 - h 7 -71 - t 3 7 -72 - t 3 J -73 - t 2 ^ 1 - h 

1 *4 *3 *2 1 £2 

"£4 r dt 3 r dt 2 f dti -1 r dt 2 f dt\ 



J —70 - £3 i — 71 - £3 •/ -73 — to J 1 — £1 J J to J 1 



— 72 —t^J —7l — £3 J —73 — t 2 J \ — t\> J to J 1 — £l 

' d ' ' 

1*3*2 1 *3 *2 

<i£ 3 /■ d£ 2 / ^1 /" ^3 /" ^2 /" c^i 



' 

d£ 4 / d£ 3 / dt 2 f dti f dt 4 f dt 3 r dt 2 f dti 



-7i - *3 J -72 - *3 J h J 1 - h J -72 -t 3 J -71 - h J h J 1 - ti 



In order to get the expression under the sign of the limit in Eq. (f2*K|) one applies the 
definition (JZj) for the multiple polylogarithms in Eq. (|2fijl . Using the same definition for 
the final multidimensional integrals in (|2*Hj) and making the change 71 — > a 2 , 72 — > 0:3 
one finally arrives at the result for the case et\ — 1 and 0:4 = 0, 

(1, a 2 , a 3 , 0) = -Lion f-a 3 , — , — — ) - Lz 2 ,i 1 f-a 2 , — , — — 

V a 3 a; 2 / V a 2 03 

- ln(a 3 + l)Lz 2 1 f-a 2 , — — ) - ln(a 2 + l)Li 2 a (-a 3 , — - 

V a 2 J V a 3 , 

-ln(a 2 + l)hi(a 3 + l)C(2). (29) 



(28) 



2. a-i = 1, a 2 = a 3 = 

For these values of the parameters one has an integral of the very simple form 

ln(l - y) In y 



L_ ++ (l,0,0,a 4 ) = / dy- 



a 4 + y 

u 

After a change of variable y — > 1 — t one gets 



? ta h,«ln'(l-«) = _> h'(l-«,) }j, = 

«4 + 1 — £ J «4 + 1 — £1 J t 2 

ti 

£ 2 •> 1 «4 + 1 - k ~ J t 2 J a 4 + 1 - *i J 1 - * 3 J 1 - t 4 ' 





Applying the definition (0) one obtains 



L_++(l, 0, 0, « 4 ) = -2Li 1>li2 f 1, a 4 + 1, ) . (31) 
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3. a% — 1, a 2 = (and a 4 = — 1 J 

We shall again find the limit of the rhs of (|26J) for t m = 1 and 71 — > 0. The first 
and the third terms are equal to because of the limit lim 71 ^o m (7i + 1) = 0. The 
other terms transform into 

lim .14,1,1 (-73,—, — ) = ~Li h 2 (-73, — -J , 
7l ^o ^ 73 71 / \ 73/ 

y r . ( 72 71 ~tm\ T- ( 72 -l\ 

hm L«i,i,i,i -73, — , — , = -Li 1X2 -73, — , — , 

V 73 72 7l / V 73 72 / 

r t- ( 7i 72 -t m \ T . ( 72 -1\ 

hm Lii,!,!,! -73, — , — , = -Wiai -73, — , — • 

7l ^o y 73 7i 72 / V 73 72 / 

Finally we write 

(1, 0, at 3 , a 4 ) = -Li 112 (-a 4 , — , — ) - Lii (2 ,i f-a 4 , — , — 

V «4 «3 / V «4 «3 

- ln(a 3 + 1)^1,2 f-«4,-—V (32) 

V Q!4/ 

For the special case 0-4 = —1 one gets 

L-++ (1, 0, a 3 , -1) = -Lii ,1 , 2 f 1, -a 3 , — ) - L*i )2 ,i ( 1, -a 3 , — 

V 013 / V 0:3 

-ln(a 3 + l)C(3). (33) 



4. «2 = «3 = ('and a 4 = — 1 ^ 

In this case one proceeds along the following lines: 

L. ++ (a 1 ,0,0, a4 )=/d i / ln(a '- i/)ln2!/ rS- /^-l+Ol^l-i) 
J 0:4 + V J 

* 

-U h2(1 - fi) / *» =-/*, ln2(1 - f '> 

7 — «4 — 1 + ti J «i — 1 + to J ~ a 4 — 1 + tl 

-ai+2 

' dt 2 /* ln 2 (l-ti) 




, ln"( 1 -t x ) . f , in (1 - 

atx Inai / dti 

aii — 1 + h J — «4 — 1 + ti J — «4 — 1 + ii 
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/ dh J ^ ?-*L_ ?^-21na 1 Li 3 (- 

J a 1 -l + t 2 J -Oi A - 1 + ti J 1 - t 3 J 1 - U V 





Using the definition (JJJ) we arrive at the result 

L_ ++ (ai, 0, 0, a 4 ) = 2Lzi i x i f 1, a 4 + 1, - — rr, — ) ~~ 2 ln a i Li 3 ( — _ ) • ( 35 ) 
For the case a 4 = — 1 one obtains 

0, 0, -1) = -2Li 1>2>1 (l, 1 - a u ) - 2 In «! C(3) . (36) 

V 1 — CCi/ 



5. a 2 = ('and a 4 = —1) 

For this integral we change the integration variable y — > 1 — t, 

1 ln(ai -y) ln7/lnQ 3 + y) = /■ ^ln(l - i) ln(ai - 1 + t) ln(a 3 + 1 - t) 



a 4 + ?/ j a 4 + 1 — i 

i 



o 



^ ln(l-t) ln( 7l + t)ln( 72 -t) 

7 3 -t ' 1 ) 



One notes that the last integral is an analogue of the integral in Eq. (J26j) . The 
calculation proceeds in a similar way, 

/ 72 71 1 \ /" / ^3 ? <& 2 / 

- ^1,1,1,1 7s, — , , = / — / / / z 

V 73 72 7i/ 7 7i + h I 72 - *3 7 7s - h J 1 - *i 

dt$ f dt s f ln(l — 1 2 ) 7 dt± f ln(l — 1 2 ) f dt s 



dt 2 ■ — = - / / dt 2 



{] 7i + U J 72 - *3 J 73-^2 J 7i + U J ' 7s - h 7 72 - h 

/rft 2 ln(1 - t2) - /-*L_ } 4 ^ln(l-t 2 )ln( 72 -t 2 ) 

() 7l+^4 7 73-^2 J ll+t 4 J 73-^2 

ln(l - t 2 ) ln( 72 - i 2 ) r dt 4 



Y'jl) - f dt 2 - — (72 - — / 



73 - *2 / 7l + *4 

r(i) - in (7l + 1) / dt 2 Hl ~ h) 111(72 ~ h) + } dt Mi-h)Mii + h)Mii-h) 

7 73 -^2 7 73 -^2 
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Y'(l) - Y"{1) + J dt 2 



ln(l - t 2 ) ln( 7l + t 2 ) ln( 72 - t 2 ) 



73 - h 



where we have introduced the notation 



/ 7! + i 4 7 73-^2 



y"( t „,)=lnfa + U/^ " l(1 " f2 " n ! 72 " f2) . (39) 

7 73-^2 

The last term in ()38|) is the required integral. The expansion of the integral Y'{t m ) 
in terms of multiple polylogarithms is similar to the evaluation of I'(t m ) in Eq. (|22|) . 
The result of the calculation is 

Y\t m ) = Lii 1 1 1 f 7 3, --, --, -) + ln( 7 2 - t m )Lix Xl L, --, ■ (40) 

V 73 7l 72 / \ 73 7i / 

For the calculation of Y"{t m ) one can make use of (fT8j) . Finally using Eqs. (}38j) and 
one arrives at the result 

1 ln(l-t)ln(7i + t)ln(72-t) _. / 7i 72 1 \ 

«* 7 = -^1,1,1,1 73, , , — 

73 - 1 \ 73 71 72 / 

-Lii,i,i,i (73, — , -— , — J - hi(7i + 1)^1,1,1 (73, — , — J (41) 
V 7s 72 71/ V 7s 72/ 

- ln(7 2 - l)^i,i,i (73, - — , -— ) - hi(7i + 1) ln( 72 - l)Li l x (73, — ) . 

V 73 7i/ V 73/ 

To obtain the formula for the L function with a 2 = we must only change 71, 72, and 
73 to cti — 1 , 0:3 + 1, and «4 + l according to Eq. (JHTj) : 

]n(ai-y)\ny\u(a 3 + y) 







L_ ++ («i, 0, a 3 , a 4 ) = / dy 



a A + y 



1-ai l + a 3 1 \ / 1 + a 3 1 - «i 1 

-^1,1,1,1 1 + a 4 , — , , — - ^1,1,1,1 l + a4, 



1 + a 4 1 — «i 1 + « 3 / V 1 + a 4 1 + a 3 1 — «i 

motiLii^x 1 + a 4 , , I — ma^Liin 1 + a 4 



1 + CK4 1 + «3 / ' ' V 1 + a 4 1 — «i 

lnailna 3 L2u ( 1 + a 4 , — - — ) . 

V 1 + a 4 / 
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For the case = — 1 we calculate the limit of the rhs of (|42j) for a± — ► — 1 and obtain 

1 + a 3 1 

1 — Qti 1 + «3 



L_ ++ (ai, 0, a 3 , -1) = Li 2 ,i,i 1 - «i, 



+Li 2 ,i,i (1 + 013, r~ — ~ » — ) + ln«iLi 2 i (l + a 3 , 



l+«3 1-fti/ V 1 + ^3 



+ lna 3 Lz 2i i ( 1 — ax, ) + ln«i lna 3 ((2) . (43) 

1 — Oil. 



IV. TRANSFORMATION OF L ++ , TO MULTIPLE POLYLOGARITHMS 



In this section we will show that all our functions can be expressed in terms 

of multiple polylogarithms. 



A. General case for the L +++ function 

We now proceed with the transformation of the triple index function 



i 



T , \ f i ln («i + v) ln ("2 + y) ln(a 3 + y) 

L +++ (a 1 ,a 2 ,a 3 ,a 4 )= dy ■ . (44) 

J otA + y 

o 

After changing the integration variable y = —a%t we obtain 

ln(ai ~~ a it) hi(a 2 — ait) ln(a 3 — ait) r lnai ln(a 2 — ait) ln(a 3 + ait) 



9L-t J SL -t 

ln(l - t)[lnai + ln(g - t)][lnai + ln(g - t)] = 

Q l 

} ln(a 2 + y) ln(a 3 + y) 2 V" 1 ln(l - t) 

lnai / dy ■ In ai / dt— — 

J a 4 + v J 21 — t 

o y o Q i 

, "f. Ml-tlMg-*) . "V /ai ln(l - t) ln(g - t) 
-lnai y dt — — lnai y rft — — 

ln(l - t) ln(sa _ t) Wsa _ t) 
dt— --^ '—^ . (45) 

OCA £ V I 



The first integral on the rhs of (J43j) has been calculated in Eq. ()19j) . For the second 
integral one makes use of the formula (fTTj) (the only change is a 2 — > ai). For the 
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evaluation of the third and fourth integrals one uses Eq. (|18|). We are left with the 
most complicated fifth integral. Let us consider an integral of the type 



This integral is an analogue of the integral in Eq. (j2HJ). The calculation proceeds in 
a similar way. One obtains the result 

^ >(l-t)ln( 7l -t)ln( 72 -t) / t m \ 

ere = - ln(7i - t m ) 111(72 - t m )lA lf i 73, — 

73 - t V 73 / 

- ln(7 2 - t m )Li llx I 73, — , — J - hi(7i - *m)^'i,l,i ( 73, — , — J (47) 
V 7s 7i/ V 7s 72/ 

j- ( 12 li t m \ ( 71 72 t m \ 

-^1,1,1,1 73, — , — , — - ^1,1,1,1 73, — , — , — • 
\ 73 72 7i/ V 73 7i 72 / 

Taking into account everything mentioned above for Eq. (|45jl we arrive at the final 
result for the L +++ function, 

/ a 4 a 2 a 3 1 
L +++ («i, a 2 , «3, «4j = ^1,1,1,1 — , — , — , 

VCUl «4 «2 «3 
/ «4 «3 C*2 1 \ / «4 «3 1 

+^1,1,1,1 — , — , — , + maiLzi^i — , — , 

\ai OL\ «3 «2/ V«2 «4 «3 

+ ln(l + 012)^1,1,1 (—> — , — M + ln (! + 03)^1,1,1 (—, — , — - 

\ai «4 0:3/ VOJi «4 O^, 

+ In ai In a 2 £*i 1 f — , — —) + In ati ln(l + a 3 )Lzi 1 f — , — M (48) 
\a 3 a 4 / \«2 «4/ 

+ ln(l + a 2 ) ln(l + a 3 )Li\ 1 f — , ^ + In at\ In a 2 In a 3 In f ^ 4 \ _ 

' \ai QI4/ V 04 / 

For this equation the conditions (jHJ) are assumed. We emphasize that the arguments 
of L +++ functions occuring in the actual calculation of the massive scalar one-loop 
integrals are not of the most general type as assumed in the derivation of (|48j). We 
have nevertheless included a discussion of the general case because Eq. ()48j) may be 
useful in other applications. In the results for the massive scalar integrals one has 
only the special cases where a\ = a 2 or a,\ = a 3 as well as the cases a\ — and/or 
a.2 = and/or 0-3 = and/or 04 = { — 1,0}. If some en's coincide with each other 
Eq. (|48J) becomes simpler. In this case one can also make use of symmetry properties 
to obtain simpler relations between the L +++ functions and multiple polylogarithms. 
For the cases 07 = and/or a 2 = and/or a 3 = and/or a 4 = {—1, 0} the general 
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formula (}4*Hj) is no longer valid and these cases must be studied separately. 



B. Special cases for the L +++ function 

In the Laurent series expansion of the massive scalar one-loop integrals the fol- 
lowing special cases for the a^ are present. 



1. ot\ = a<i or ai = 03 

As it was stated in Sec. II the L +++ function is symmetric under the permutations 
o.j <-> (Xj. Therefore, it suffices to consider the case Oi = o 2 . 
We must evaluate the integral 

t i \ } j ln2 ( a i + y) ln ( a 3 + y) , Atis 

L +++ (ai,ai,a 3 ,Q!4) = / dy ■ . (49) 

o 

This integral can be expressed in different ways. First of all one can directly use 
Eq. (|38j) replacing o 2 by a.\ . The second possibility is to use symmetry properties. 
One takes into account the rhs of Eq. (|48|) and notes that the part with multiple 
polylogarithms of weight four is symmetric under the exchange o 2 *-> 03. It allows 
one to reduce the number of the multiple polylogarithms from two to one. First we 
apply Eq. (}4*%|) for the case o 2 = 03 replacing 03 by <x 2 . Second we change a\ — > 03 
and a.2 — *■ oc±. After these transformations one obtains the following result: 



± 

L +++ (a t , 01,03,04) = J dy 



ln 2 (oi + y) ln(o 3 + y) 
a 4 + y 



+2^1! ! ! — , — , 1, + In o 3 Li\ \ i — , — , 

V03 O4 Oi\) \OL\ 04 Oi 

+2 ln(l + a x )Li x hl (—, —, — -) + In o 3 [ln(oi + 1) + In o x ] Li 1A (—, — - ) (50) 
\a 3 a 4 a,\j V01 04/ 

+ ln 2 (l + a\)Lii 1 ( — , + In 2 a\ In 03 In ^ 



03 04/ v 04 

There is also the third possibility to express L +++ (a\, 01,03, 04) in terms of multiple 
polylogarithms: 

^ In (01 + y) ln(o 3 + y) y ^-_ ai t f ^ln (01 - ait) ln(o 3 - ait) _ 



1/01 

dt 



In 2 01 + 2 In 01 ln(l -t)+ ln 2 (l - t) In a x + ln(^ - t) 



04 £ 

«1 
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+ lnV l^-r + l^a, /^^-21nV ^ dt^^- 



«4 + y J «4 + y J — — t 

o y o y o "i 

_ lnai r d( ^_ 21nQ1 " i /v (i - f)in( »- i) 



«4 _ t j «4 _ t 



ln*(l -t) ln(^ -t) 
dt J_] ai (51) 



The first term can be integrated immediately. For the second and third term one uses 
Eq. (JTfij) and Eq. ([17)1 . respectively. The integral of the fourth term can be rewritten 



as 

-l/ctl , 2/- .s ~l/ a l *1 ,. *2 



"i «l 

The fifth term is calculable with Eq. ()14j) . To integrate the last term one first evaluates 
the following integral: 

K ln 2 (l-t)ln(A-t) = _ V m rfti ln 2 (l - h) \ ) , J \ dt 



fo-t J fc-h ]J i A-*2 

U [tm pi — 1 J 

>" rft 2 ln^l-tO , , fl ln 2 (l-t) 

2^i,i,i,i (^'f >|^) +21n(A-0^i,i,i • 

Then to calculate the last term of Eq. (JST]) one only has to change /?i,/3 2 , and t m 
by the corresponding combinations of a^. Finally we arrive at the result for the 
ai, a 3 , a 4 ) function, 

/ «4 tt 3 1 \ - , / «4 1 

L +++ (ai, a 1( a 3 , a 4 ) = -2Lzi 1 1 a 1, — , — , 2 ln(a 3 + 1 1 1, — , 

V «i «4 03/ V «i a 4 

/ «4 a 3 1 \ - . / «4 1 

+21naiL^i 1 1 — , — , + 21n«i ln(a 3 + l)Lt 1 1 — , 

\Oii «4 a 3 / \Oti «4 



1 2 r- ( a i 1\ ,2 1 1 /«4+l\ 

+ in a.\L%\ 1 — , + m a.\ In a 3 In . (54) 

' Va 3 0L4J V «4 / 



This is the third possibility to express L +++ (a\, a%, a 3 , 0:4) function in terms of multi- 
ple polylogaritms. Each of the Eqs. (jSUJl and (p>4"j) contains only one multiple polylog- 
arithm of weight four and they are both equally acceptable from this point of view. 



17 



One has a free choice to apply any of these equations for the required L functions. The 
situation with the a%, a 3 , a 4 ) function is an example of the statement that 

the expansion of the L functions in terms of multiple polylogarithms is not unique. 



2. «i = (or a 2 = or a 3 = 0) 

For this integral we change the integration variable y — > 1 — t, 



i 



.(0, a 2 , a 3 , a 4 ) = J dy 



\\iy\\i(a 2 + y) ln(a 3 + y) 

a 4 + y 



o 



1 ^ ln(l - t) ln(q 2 + 1 - t) ln(a 3 + 1 - t) 



and using Eq. (|^8|) we arrive at the result 

1 + a 2 1 + « 3 1 



L +++ (0, a 2 , a 3 , a 4 ) = 1 + a 4 , 



1 + a 4 1 + a 2 1 + «3 - 



/ n 1 + a 3 1 + a 2 1 \ . r . A 1 + a 3 I 

-^1,1,1,1 1 + "4, ' , — , — - lna 2 ^i,i,i 1 + ^4, 



1 + a 4 1 + a 3 1 + a 2 / V 1 + a 4 1 + a 3 / 

In a 3 Li 1AA (l + a 4 , - + " 2 , — - — ) - In a 2 In a^L^i f 1 + a 4 , — - — ) . (56) 
V 1 + a 4 1 + a 2 / V 1 + a 4 / 



5. ai = a 2 = 

To calculate this integral we again change the integration variable y — > 1 — t, 

i „ 

In y\n(a 3 + y) 



L +++ (0,0,a 3 ,a 4 ) = / dy- 



a A + y 



I h'd-QMc + i-t) 

7 04 + 1 - t 

o 

For the last integral we use Eq. (J53j) . An additional simplification can be done if one 
notes that 



Li lxl [l,a4 + l,—!— ) =-Lig(-— ). (58) 
a 4 + l/ V a 4 / 



Finally one has 



ct 3 + 1 1 \ / 1 . 

L +++ (0, 0, a 3 , a 4 ) = 2Li x x 1 1 1, a 4 + 1, — , — - 2 lna 3 Li 3 . (59) 

a 4 + 1 a 3 + 1/ V a 4 / 
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4- «l = «2 = 0, «4 = —1 (or (X 2 = (X 3 = 0, «4 = —1) 

In this case one should calculate the limit of the rhs of ()59|) for t m = 1 and 
0:4 — > —1. After this procedure one obtains 

L +++ (0, 0, a 3 , -1) = -2Lzi,2,i fl, a 3 + 1, r ) -21na 3 C(3). (60) 

For the case a-2 = a 3 = and 0:4 = —1 one can use the same formula. The only 
change is 0:3 — > a\. 



5. cui = 0, «4 = — 1 

To obtain the solution for these values of the oti we must find the limit of the rhs 
of ()56|) for 0-4 — > —1. After taking the limit one arrives at the result 

/ 1 + ^3 1 

L +++ (0, a 2 , a 3 , -1) = +Li 2X1 [l + a 2 , 



1 + a 2 1 + a 3 
+£22,1,1 (1 + 0:3, : + ^ , z ) + lna 2 Li 2>1 fl + a 3 



1 + a 3 l + a 2 J ' ' V l + a 3/ 

+ lna 3 Li 2 ,i ( 1 + a 2 , ) + lna 2 lna 3 ((2). (61) 

V l + a 2 / 



V. TRANSFORMATION OF L + TO MULTIPLE POLYLOGARITHMS 

In this section we will show that all our L + functions can be expressed in terms 
of multiple polylogarithms. 



A. General case for the L + function 

Here we derive the general formula for the single index L + function Eq. (J2J), 

L + («i, a 2 , a 3 , a 4 ) = / dy ^—^—^-Li 2 (a 2 + a 3 y). (62) 
J otA + y 


After changing the integration variable y — > (t — a 2 )/a 3 one gets 



■ # Hai + ^) = I rff -'"°3 + ln( aiQ 3- Q2 + t ) L . 
a 3 0:4 H 22 j a 3 tt4 — a 2 + t 

«2 a 3 «2 
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The integration interval can be split into two pieces, [0:2, 0] and [0, a 2 + a&j. One can 
then write L + as a sum of four terms, 



CT2+03 02 ^ , ( Q 2 



— -I / -/[-^L i2(t)+ /-/ <£££W (64, 





where we have introduced the notation 



1+t J J \ j + t 



a = a\as — a 2} 7 = 0304 — a^- (65) 
Looking at Eq. ()64|) it is clear that there are only two different types of integrals to 



Li a (t) and / dt ^ + j Li 2 (t). (66) 



be dealt with, 



7 +t 7 7 +t 

u U 

The upper limits are t m = a 2 + 03 or £ m = o: 2 . The first integral can be evaluated 
analytically in terms of standard logarithms and classical polylogarithms up to Li 3 . 
However, the same integral can also be expressed in terms of multiple polylogarithms 
via the integral representation (J7J), e.g., 

'"' * Li2 w = "[Jh r 'f d Ji'fj!h r = _ Llu f_ 7 , Zhl) . (67) 



7 + t J ^ + ti J t 2 J I -h ' V 7 

' 000 



We now deal with the second integral in (J66J). Consider the following multiple poly- 
logarithm of weight four: 

T . ( a t m \ ? dt 2 J dti y . / , \ y dti f J dt 2 

Li2,i,i -7, -, — = / / -Li 2 (ti) = / — — Li 2 (ti) / — — - = 

\ 7 -a J J Q -a - t 2 J -7 - h 7 1 t a + 2 

*? dt *? dt 

/ -^-Li 2 (ti) ln(« + t m )- / _J-Li 2 (ti) ln(a + h). (68) 

' 
In the first step we have used the usual trick to change the order of integration. As 
already noted before [see Eq. (|57jl] the first term on the second line can be expressed 
through a multiple polylogarithm of weight three. Thus one has 

"J dt Ma±Au 2 {t) = -Li 2jl>1 f-7, °, *") - Li 2>1 f-7, In(« + O- (69) 



Finally, substituting Eqs. (}6T|) and (JBT?)) into Eq. we arrive at the desired relation 

a 2 — a±a^ a 2 



L + («i, a 2 , 0:3, 04) = Li 2 ^\ ( a 2 — «3«4 



«2 ~ O3O4 ' «2 — ttl03 



T ■ ( "2 - «ia 3 «2 + «3 \ . , r • / "2 \ /- n v 

-Li 211 a 2 - a 3 a 4 , , + Ina^jy a 2 - 03^4, (70) 

V «2 — «3«4 «2 — «1«3/ V «2— «3«4/ 

- Infai + 1)L? 2 1 a 2 - 0304, • 

V a 2 — a^a^J 
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We should note that, similar to Eq. (|27j). the conditions (J5J) are assumed for the 
variables a^. Also, one cannot directly use Eq. (fTUj) if o 2 — a 3 a 4 = or o 2 — 0:10:3 = 0. 
However, in the results for the massive scalar integrals precisely these special cases 
appear, as well as the cases where ai = and/or o 2 = and/or 03 = and/or 
04 = { — 1,0}. In such cases the general formula (fTT)|) is no longer valid and these 
cases must be studied separately. 



B. Special cases for the L + function 

In the Laurent series expansion of the massive scalar one-loop integrals the fol- 
lowing special cases appear for the arguments of the L + functions: 



1. ot-2, — 0304 = (or 0(2 — 010:3 = 0) 

In this case one must find the limit of the rhs of Eq. (fTOj) for o 2 — > «3«4- First we 
rewrite the rhs of Eq. (|7U|) in terms of multidimensional integrals via the definition (JJJ). 
Second we replace o 2 by 0304. We finally again use the definition ((7J) to obtain the 
result 

L+(oi, a 3 a 4 , 03, "4) = -Li 3 ,i ( o 3 (o 4 - «i) 



04 — Oi\ 



+Li 3 1 fa 3 (a 4 - Oi), ° 4 + — ) - lnoiLi 3 (a3a4) + ln(oi + l)Li 3 (03(04 + 1)) . (71) 

V O4 — OL\) 

When o 2 — Oi0 3 = one must find the limit of the rhs of Eq. (|70jl for o 2 — > O1O3. 
We again rewrite the rhs of Eq. ()70|) in terms of multidimensional integrals. We then 
replace o 2 by 01O3 and use the definition (JJJ). We arrive at the result 



L + (oi, oio 3 , o 3 , o 4 ) = —Li 2: 2 [o!3(ai — o 4 ), 
+Li 2j2 ( o 3 (oi - o 4 ), ai + 1 ) + lnoiLz 2j i fo 3 (oi - o 4 ), — — — ) (72) 



Oil ~ ^4 

- ln(oi + l)Li 2 ,i ( ol 3 (o£i - o 4 ), 



Oil 


Oi 




o 4 




Oil 




Oi\ 




o 4 


a l 


+ 


1 N 



Oi — O4 



2. 01 = 

Unfortunately in this case one cannot use Eq. ()70|) for 01 = because one is im- 
mediately faced with the problem of a logarithmic infinity. One must find another al- 
gorithm to express the L + (0, o 2 , o 3 , o 4 ) function in terms of multiple polylogarithms. 
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After changing the integration variable y — > 1 — t one gets 



/ dy ■ — L12O2 + a 3 y) = [ dt — -Li 2 (o 2 + a 3 - o 3 t) 

J «4 + y J «4 + 1 — t 



1 1/1 1 \ *i 



; _ ; ln(l-t x ) > ^ ln(l - a 2 - o 3 + a 3 t 2 ) 



04 + 1 - ti 7 sa + 1 _ t 



«2/o:3+l 



«3 



2 



1 0:2/0:3+1 

, ln(l - a 2 - a 3 + a 3 t 2 ) f , . ln(l - ti) r/ u . / 1 
cft 2 — / dti — — - Li 2 (a 2 )^i,i a 4 + 1 



% + 1 - h J —a A + 1 - t x ^ V* ' a 4 + 1 

The last integral is an analogue of I'(t m ) in Eq. (|22|). First one notes that 

V M ^7^ = "^M («4 + 1, -XT) • (74) 
o 4 + 1 — ti V o 4 + 1 / 





Then one considers the following chain of transformations: 

1 *2 

/ / 7 ^1 1 I ola + 1 

7 1 - o 2 - o 3 + o 3 t 2 / ^ + 1 - t! 1,1 V 



, OtA + 1 

^ J ^ a 3 

J g + l-ti V an + lJJ I - a 2 - a 3 + a 3 t 2 



~ ln(l - a a ) / , 24 ltl (o 4 + 1, (75) 

a 3 J g + 1 - tx V 0:4 + 1/ 



1 /" ln(l - a 2 - a 3 + a 3 ii) . / t 4 

«*1 m 7~, T ^1,1 « 4 + 1, 



a 3 y ^ + 1 - *i v «4 + 1 



«3 



Using Eq. (|74|) we see that the last integral is exactly the integral required in Eq. (|73|) . 
The initial integral of Eq. (J7H|) and the first integral of the rhs of Eq. (J7H|) can be 
expressed in terms of multiple polylogarithms due to the definition ( |7j) . Finally for 

the L + (0, a 2 , a 3 , a 4 ) function we obtain 

T (n \ t ■ f a 2 + a 3 a 2 + a 3 -l a 3 \ 

L+(0, a 2 , a 3 , a 4 ) = "4 + 1, — 7 — r, : , : - (76) 

\ 03(04 + 1) a 2 + a 3 a 2 + a 3 — 1 ) 

+ ln(l - a 2 )Li 1)lt i ( o 4 + 1, — ^r-^-^rr-, — — — J - Li 2 (o; 2 )ivi 1| i fo 4 + 1 



3 (o 4 + 1) ' «2 + O3 / 'V 'a4 + l 
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For these values of the ai one uses Eq. f!76|) to calculate the limit of the rhs for 
a 4 — > —1. One arrives at the result 

T fn 1\ T- /«2+«3 «2 + «3-l «3 \ 
L+(0, « 2 , « 3 , -1) = -^2,1,1 > j > j 7 ( 77 ) 

V a 3 a 2 + a 3 a 2 + a 3 — 1 / 

- ln(l - a 2 )Lz 2il C^±^ f — ^— ) + Li 2 (a 2 )C(2). 
V a 3 a 2 + a 3 / 



4- «l = 0, «2 + «3 = 1 f ^nd « 4 = —1 ) 

If one takes a look at Eq. (J76j) one realizes that there is a problem if a 2 + a 4 = 1. 
To express the L + function for this configuration of the oti the limit of the rhs of (|76j) 
for a 2 — > 1 — a 3 must be found. The result is 

L+(0, 1 - a 3 , a 3 , a 4 ) = -£21,1,2 ( a 4 + 1, — r~~~7T > a 3 ) (78) 

V a 3 ( a 4 + 1) / 

+ Yn.a 3 Li l lx (a 4 + 1, — , 1 ■ a 3 J - Li2(l - "3)^1,1 f«4 + 1, — ^— 7) ■ 
\ a 3 (a 4 + 1) / V 0:4+1/ 

For the case ot\ = 0, a 2 + a 3 = 1, and a 4 = — 1 one must find in addition the limit 
for a 4 — > —1. One arrives at the result 

L+(0, 1 - a 3 , a 3 , -1) = Li 2 , 2 (—, a 3 ) - mo 3 Li 2 ,i f — , o 3 ) + ((2)Li 2 (l - o 3 ). (79) 

\0 3 / \0 3 / 



5. (Hi = 0; a 2 = —oi 3 

To obtain the result for this case one must calculate the limit of the rhs of (fTE|) 
for a 3 — ► — a 2 . After taking the limit one has 

L+(0, o 2 , -a 2 , a 4 ) = - £01,1,2 f — ^-7, -a 4 , — — ) 

\a 2 — 1 a 4 / 

77) • (80) 
a 4 / 







f-04, -) 


+ Li 2 (a 2 )Li 2 ^ 


V a 4 / 
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6. «i = 0, o 2 = 

For this case one can directly use Eq. (|7S|). 

(1 O3 — 1 O3 \ 
"4 + 1, — ' ' T ■ ( 81 ) 
O4 + 1 O3 « 3 — 1 / 

But there is also another very simple possibility. We first change the integration 
variable y — > t/a 3 , 

j dy J^U 2 (a 3 y) = /*^Li 2 ( t ) = = 

J at4 + y J 0:30:4 + t J 03O4 + ti J fa 

«3 

/ / dt, = rdfar dt, rdfarjt^ 



OL-iOiA + t\ J to J —Ot-iOLA + £1 7 £3 J 1 — £4 

" d ,i4 i o d o 

Now using the definition ((7j) we obtain the result 

L+(0, 0, o 3 , o 4 ) = Li 2 , 2 f-a 3 tt4, —) ■ (83) 

V o 4 / 

The reader has a free choice to use either formula (|81|) or ()83|1. Both equations contain 
multiple polylogarithm of weight four. The depth of the multiple polylogarithm in 
Eq. is two against four in Eq. ()81|). For o 4 = —1 Eq. can be directly used. 
However, in the case of Eq. (JHTJ) one must first calculate the limit for o 4 — > — 1. 



7. Oi = 0, o 2 = 1 

Unfortunately, in this case one cannot use Eq. (fTHJ) because of the term ln(l — 02). 
To express this L + function in terms of multiple polylogarithms we first make use of 
a standard relation between dilogs with arguments x and 1 — x for the function Li 2 
under the sign of the integral: 

f In y f In y 

/ dy ; — Li 2 (l + a 3 y) = / dy ■ — [C(2) - ln(-o 3 y) ln(l + a 3 y) - Li 2 (-a 3 y)] = 

J 04 + y J 04 t" y 



\ny } hxy . f hiy[hx(-a 3 ) + hiy] ln(l + a 3 y) _ 



C(2) f dy—^- f dy — U 2 (-a 3 y) - f dy 

J OIa +v J 04 + y J 



o 4 + V J otA + y ~ J q o 4 + y 

C(2)Li 2 (--) - Li lxl>1 U + 1, — |— , -^t) (84) 

V 04/ V 04 + 1 03 03 + I/ 

lnyln(l + o 3 y) /■ ln 2 ?/ln(l + a 3 y) 



H- a3 )[dy lavln[ \ +a3V> - [dy 

J O4 + W J 



04 + y J 04 + y 

24 



where the £#1,1,1,1 function was obtained with the help of Eq. (|81j). To obtain the last 
integral in Eq. (J8%j) one proceeds as follows: 



^ ln 2 yln(l + a 3 y) „_^_ t > ^ln 2 (l - *) fn(l + a 3 - a 3 t) 



«4 + y 7 a 4 + 1 — t 







a 4 + 1 - £1 •/ H-CK3- "3^2 / + 1 - *2 / 04 + 1 - h 

u 1 " 3 

? _dh__ r dti r dt 3 r dt 4 
J — + 1 -t 2 J a 4 + 1 -ti 7 1 -t 3 J 1 -t 4 ~ 

Q 3 

or- fi , 1 "3 + 1 «3 \ 

2L«i,i,i,i 1, a 4 + 1, — - ( — -, — ■ 

V a 3 (a 4 + l) a 3 + lj 

Similarly one can evaluate the remaining integral 

1 lnyln(l + o 3 y) / a 3 + 1 a 3 \ , . 

"3/ ; = --C«i,i,i a 4 + l, — 7 j TV) — r ■ (86) 

a 4 + y \ a 3 (a 4 + 1) a 3 + 1 / 

Now combining Eqs. (JBU), and one arrives at the result 

"3 + 1 "3 \ 



L+(0, 1, a 3 , 04) = -2Lii,i,i,i 1, a 4 + 1, 



03(04 + 1) a 3 + 1 



/ / , 1 1 a 3 + 1 a 3 , 

-L«i,i, M 04 + I, — , , — 87 

a 4 + 1 a 3 03 + I/ 



+ ln(-o 3 )^i,i,i ( "4 + 1, y~~ ~~ ~TT ) -^-r ] + C(2)Li 2 (-— ) . 

\ o 3 (o 4 + l) a 3 + 1 j V a 4 / 



OL\ = 0, «2 = — O3 = 1 

For these values of the a, we must find the limit of the rhs of Eq. (JHTj) for a 3 — ► — 1. 
After taking the limit we obtain 

L+(0, 1, -1, o 4 ) = Lii 1 2 (o 4 + 1, — 77, 

V o 4 + 1 / 

+2^1,1,2 f 1, o 4 + 1, —^—) + C(2)Li 2 (-—) . (88) 
V o 4 + l/ V 04/ 



VI. CONCLUSIONS 

We have presented all the necessary relations to transform the L functions [as 
defined in Eqs. (JTJ) and (J2J)] that occur in our 0(e 2 ) results [I] for the Laurent series 
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expansion of massive scalar one-loop integrals to multiple polylogarithms. We have 
used these relations to transform our results on massive one-loop integrals involving 
L functions to corresponding results involving multiple polylogarithms. The multiple 
polylogarithms results are readily available in electronic form [TT] . 

Despite of the fact that the relations between the L functions and the multiple 
polylogarithms have been derived having the massive scalar one-loop integrals in mind 
they can also be used in a more general setting. In fact, any definite integral given 
by 

B B 

r m(ai + b\x) ln(a,2 + &2^) hi(a3 + b 3 x)dx r ln(ai + &ia;)Li2(a2 + b^x)dx 
J a 4 + 64X J a 3 + 6 3 x 

A A 

can be written in terms of multiple polylogarithms with the help of the relations 
presented in this paper. It is worthwhile to mention that all the equations presented 
in the present paper have been also checked numerically. 

We have found several examples where the representation of the L functions in 
terms of multiple polylogarithms is not unique. This reflects the fact that multiple 
polylogarithms obey quasishuffle and shuffle Hopf algebras and hence satisfy numerous 
identities as is the case for the classical polylogarithms. More information about 
identities between multiple polylogarithms can be found, e.g., in [5] and P, and 
references therein. 

For future parton model applications of our results numerical efficiency is an im- 
portant issue. We are presently writing numerical C++ codes to compare the nu- 
merical efficiency of the two representations in terms of L functions and multiple 
polylogarithms. 
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APPENDIX 

In this Appendix we consider as an example the real part of the 0(e 2 ) coefficient 

(2) 

ReZ^i of the Laurent series expansion of the massive box Di with three massive 
propagators. Using the rules written down in the main text of this paper we have 
expressed the corresponding results of [!_ involving L functions in terms of multiple 
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polylogarithms. The L function structure of Re D\ in [TJ is sufficiently rich to provide 
an illustration of the corresponding complexity in terms of multiple polylogarithms 
when transforming to the latter representation. We mention that all multiple polylog- 
arithms up to weight three have been reexpressed in terms of classical polylogarithms. 
We then used automatic program codes to simplify the classical polylogarithms as 
much as possible, as was also done in [T]. 

We use the notation and the conventions of PQ. In brief, we use the Mandelstam- 
type variables 



s=(p 1 +p 2 ) 2 , t = T-m 2 = (pi-p 3 ) 2 -m 2 , 



u 



U -m 2 = (p 2 -p 3 ) 2 -m 2 (Al) 



for the 2 — > 2 partonic process a(p±) + b(p 2 ) — ► Q(ps) + Q(Vi) with p\ = pi = and 



m 



One finds 
ReD? - 



''. We also introduce the abbreviations 



z 3 = (s + 2t + s(3)/2, 
z 5 = (2m 2 + t + tp)/2, 



z 4 = (s + 2t-s{3)/2, 
z 6 = (2m 2 + t-tp)/2, 



(A2) 



l s = In—, 

m z 

lp = In j3, 



In— 5, It 

m z 



In 



z?, 



I 



m c 



?A 



In- 



In 



—T 

9"' 
m z 

m 2 



L = In x, 



1 



192 



{-109Z* + 240^ + 321% + 264/ 2 / 2 - 200/ T / 3 - 177/*- (A3) 



9Ql 2 T l x l z3 + l92l T l 2 J z3 + 12l% 3 + 9Ql T lJ 2 z3 - 32l x l 3 z3 - A80l 2 T lJ z4 + 
24/ T / 2 / 24 + 180/ 3 / 24 - 144/ 2 / 23 / 24 + 480/ r U 2 4 - 336/ 2 / 2 4 - 96l x l z3 l 2 z4 + 
320/J 3 4 - 1684 + 480#U/3 - A80l T l 2 J p - 4:01% - \92l T l x l z3 l p + 
336/ 2 / z3 //3 + 96lJ 2 Jp - 3Ml T lJ z dp + 33Ql%dp + I92l x l z3 l z4 lp + 
961JU/3 + 192l 3 J p + 96l T l x l 2 + 24/2/2 _ qqij^ _ 96 ^2 + 32lj3 + 



32/^ - 32/ 3 (9/ T + 20/,, + 25/ 23 + l z4 + 131 p) - 4/ 3 (8/ T - 36l x - 9l z3 - 
43/ 24 + 94l p ) - Ql 2 s (52Z 2 - 28/^ + 15/ 2 . + 2Ql x l z3 + 4Ql x l z4 + 2ti z3 l z4 + 
32/ 2 4 - Al T (9l x + 8l z3 - 15/ 24 - tip) - QOUfs - 2ti z3 lp - hU z4 lp + 76/J + 
l t (-Ul T + 60l x - 8l z3 - 60/ 24 + 8/^)) - 24/ 2 (Al 2 T + 15/ 2 . - 8l T (2l x + 
2l z3 + 4/ 24 - 3l p ) + 4/ x (3/ 23 - 8/ 24 + 13Z/J) + 2(-4/ 2 3 + / 2 4 + 12/ 23 / 24 - 
I2l z3 l(3 - 8l z4 lp + 10/ 2 )) + 8l t (8l 3 T - 31/3 + /2 ( _ 6/z3 + 33^ + + 

Ql x (3l 2 z3 - 22/ 2 4 - 2/ 23 (5/ 24 - 9lp) + 20Up 



2 3 -22/ 2 4 - 2/ 23 (5/ 24 - 9/^ + 20/^-10/5)-' " -' 2 

3 1 n! /; i \2 io;2 ; i n; ;2 , r,;3N , c;2 



10/ 2 ) + 4(-2/ 2 3 3 + 6/ 2 3 / 24 + 



164 + 9/ 23 (/ z4 - W - 18/ 2 V/3 + 9/ 24 /^ + 21$) + 6/ 2 (3/ x - 4/ 23 + 4/ 24 - 
tip) - 3/ r (5/ 2 . - 4/,(5/ 23 - 4/ 24 ) - 8(/ 2 2 3 - 2/ 23 / 24 - 3/ 2 4 + 4/ 24 /^ + / 2 )) N 
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Al s (80/i 3 + g/ 3 + lg/ 3 _ 27l 2j z3 _ g/ 3 3 _ 165 /2 /24 _ 7 2lJ z3 l z4 + 48U 2 4 - 

24Z 23 / 2 4 - 64Z* 4 + 181% + 120l x l z3 lp + 2Al 2 z3 lp + 72l x l z4 lp + A8l z3 l z4 lp + 
72l 2 Jp - 8AIJ 2 - 24Z 23 Z 2 - 48Z 24 Z 2 + 40/J - 12Z 2 (13Z X + 2Z 23 - 6Z 24 + 
6^) + 12Z 2 (8Z T + 5l x - 2Ql z3 - 10Z 24 + 121 p) - 12Z t (5Z| + 6Z 2 - 7Z 2 3 - 
10/ 23 /,4 - 16/ 2 2 4 + Z x (14Z z3 + 27/ 24 - Alp) + lOZ^ + lU z4 lp + 2Z 2 + 
l T (-5l x - 2l z3 + Al z4 + 8lp)) + 12Z T (11Z 2 + l x (8l z3 + 7Z 24 - 4^) + 2Z 2 3 - 
6Z 2 4 - Al z3 lf3 + 8l z4 lp + 2l 2 )) } + (3/ 2 /4 - 7Z 2 /2 + l T l x + 3Z 2 /4 + 5Z X Z 23 + 
5Uz4 + 1 2 J2 + l t (2l T -l x + I0l z3 + 21 z4 - 15lp) - 8l x lp + 21, J p - 111 2 - 
Wk + h + 8l x + 5/ 23 + l z4 + 6l p )) C(2) + (Sl s + 2l x )((3) - 35C(4)/4 - 

«4(?) + ^(^) + SH.g)(-u«-^ + »e + 

8Z Z Z 24 - 24Z 2 4 + Z s (4Z t + 26Z, + 24Z 24 - Sip) + 2Al x lp + 16l z4 lp - 8Z§ - 
Al t (9l x - Al z4 + Al^)) + Li 2 _ (ll/ 2 /8 - I 2 + 2l t (l x + l z4 - 2lp)+ 

U-k + 7l x /A - l z4 + 2lp) + l x {19l x - 2Al z4 + 16/^/8) + ^U 2 {-x) x 
(-15Z 2 + Al 2 t + U-39l x + 32l z4 - 16lp) + 2Z s (10Z t + 5l x - 8lp)+ 

l t (-AAl x + 32l p )) + Li 2 (^-2U 2 + 1 (-5Z 2 + 20Z 2 - 7Z 2 - 

4ZJ 23 + 12/J.4 - 4Z 2 4 + l s {-2l t + 4(2/, + l z3 + Z 24 - Ip)) - 4/J^ + 
l t {-22l x - 8l z3 + 8/^) - 8C(2)) ) + U 2 {x) ^2Li 2 + 

2Li 2 ( "^^ - 7l 2 s /A - 3/ 2 /2 - ll/ 2 /4 + l x l z3 + 3l x l z4 + Z 2 4 - 
2l x lp + h(-Al x + 2l z3 + Alp) + l s (Al t + l x /2 - 

lz3 — lz4 — 2lp) - 6C(2)) + 

L ^-)(- 4Li2W - 4 U 2 (2)-4L i2 g) + 4L^=^) + 

/ 2 /8 - 3/ 2 /2 + / 2 /8 + 2/ x / 23 - 4/ 2 4 + l t {3l x /2 + 4/ 23 + 4/ 24 - Alp) - 

2l x lp + Al z4 lp - 2l 2 - j(18l t + 7l x + 8l z3 - m z4 + 8/^) - 4C(2)^j + 

/T\ ( I 2 9 17 \ 

Li 2 (^~J \ "8 ~ ^ f ^ s ^4^ _ _ 2lt(l x — hi) ~ lx{~glx — U)j + 

Li 2 ( ~ t{ y/ ) ) f-2Li 2 f ^ - 7/ 2 - / 2 /2 + 2/J 23 - l x l z4 - l 2 zA+ 



2m? I V V-24 



28 



l s (l t + 3l x /2 — 2Z 23 + Z 24 — lp) + 31J/3 + 2l z Jp — lp — h(lx — 4Z 23 — 2/ 24 + 
2lp) + 12C(2)) + Li 3 (^^) " 7/*) + 5Li 3 k + Li 3 x 

( 5 / s _ 6/< _ ny/2 - Li 3 (4Z t + 6Z X ) + Li 3 (3Z s /2 - 5Z, - 

7Z./2) + 4Li 3 (*) (Z s -l t - 2l x ) + Li 3 (-^p) (*. - 2l t - l x ) + 
Li 3 (^2) & ~h~ Q + Li 3 (l s + bl t - l x ) + 2Li 3 (^j l x + 

Li 3 (|) (4z s -3Z t -^ + 4Z, 4 -4Z,) + Li 3 (^^)x 

(rl ls ~ k ~ \ lx + AlzA ~ 4//3 ) + Li3 ( t(7+ff) ) (~^ s ~ 2/ * + 2/z4_ 

I - 2Z^ + Li 3 (|) (| - l t - I + 2Z 24 - 2Z^ + 2Li 3 (-x)(3Z s + 
2Z t + 2Z X - 4Z 24 + Alp) + Li 3 ^ + 4Z t - | - 2Z 24 + 21 + 
Li 3 (^Tg)) + 2it + I " + 2Z„) + Li 3 (|) (3Z S + 3Z, + 
3Z, - 4Z 24 + Alp) + Li 3 Q (^ + k+ \h - 2Z 24 + 21 + 

Li 3 (^ 2{ l~ P) ) + % + \l x ~ 4Z 24 + 4Z,) + 
Li 3 (x) (yZ s - 2Z t + h x - 6Z 24 + 6Z^ + Li 3 (- 



'T 



(-21 + 



l t - 2l x + 2Z 24 - 2lp) + 2Li 4 (x) - 4Li 4 + 4Li 4 ( j) - Li 4 
Li 4 f ^ + Li 4 f ^ + 2L 14 (^) + 3L 14 f^') + 



3Li4 (ktt^j +2 ^ 1 r^'^j- 

/ m 2 x^ 3 2T \ / ^ 

2Ll ^ {-^tjt' ) " 1,2,1 r ^ J 
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m 2 T z 5 \ ( m 2 T z & 



\ T z 5 m 2 J \ T zq m 2 J 

or- ( Z3 Zi Z A or- ( m2 i Z5 \ or- ( m 2 z 5 T 

or . (s(l-P) z 5 z 6 \ ( m 2 xz 3 sT(3 z 3 \ 

2Ll2 ^ V + 2Lt2 ' 1 ' 1 \~~aT' 7p ) ~ 

0T . ( m 2 xz 3 sTf3 zA f s(l + (3) z 6 \ 



+ z 6 (s(± + 0)_ z ± z^_ 

2z 3 >T>m 2 ) ZLl2 ^\ 2z 3 'z^m 2 

T z 6 z 5 \ ( s(l-P) z 5 T \ 



or ■ ^ ^ M^r- fi + 2 6 T 



A i I. Z A 

T Z4 t 

or,- ( t T z? > T \ or,- ( t T Z3 Za 

\T z 3 T t J \T z 3 z± t 

ft T z 4 \ nT . ft T z 4 T\ 
2L.i,i,i,i 1. T J + 2^1,1,1,1 -, y , yj + 

ft T z A z 3 \ ft z 3 T 

2Lii,i,i,i -,—,-,- + 2^1,1,1,1 — , 1, — , - 





T 


f) 


) +2^i,i,i,i 


- 1, 


^3 





!T z 4 z 3 t ) ' ' ' VZ3 T t 

( t z 3 T\ ( t z 3 T z 3 

or • f t z 3 T z 4 \ ft z 3 z 4 T 

2L%W W t ' V 7 J " 2L ' 1 ' 1 - 1 ' 1 U' T ' T 

/ f z 4 T\ T f t z 4 T\ 

2Lilw ih ¥• 7 3 - 7 ) - 2Lk ^ (h ¥• h f) + 

(t Z Z T\ fT Z Z \ 

.t'^^'TJ"^ 1 ' 1 ' 1 ' 1 1 v ^' 1 't'^J + 

(m 2 z 5 T\ (m 2 z 5 T z 5 \ 
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T . fm 2 z 6 T \ (m 2 Zq T z 5 \ 

^1,1,1,1 ) 777) - 1 ) — 7 ~~ J 777) j — o ~~ 

y ^6 i m / \ zq 1 z$ m z J 

T . (m 2 z 6 z 5 \ (m 2 z 6 z 5 T \ 

^1,1,1,1 — j — > — - ^1,1,1,1 — > — > 777' - 2 ~ 

^1,1,1,1 777, — ) 1, — 2 + ^1,1,1,1 — ) 777, 1> — - 

( 8(1-0) z 5 z 6 \ ( 8(1-13) z 5 T V 

3-^1,1,1,1 n ' — ' — 2 + ^1,1,1,1 n ' 777) — ) — 7 

\ ZZ4 zq m 2 ) \ 2z± 1 Zq m 2 

(3) z 5 z 6 T 



( s(l-(3) z 5 z 6 \ ( 8 (1 



2z4 ' z§ T ' m 2 , 
qr . ^(1-/3) z 5 z e z s \ ( s(l + /3) z 6 T N 

"5^1,1,1,1 o ) — > — > — 7 + ^1,1,1,1 7, > 777) — 7 

\ 2z 4 zq z<r> m 2 ) \ 2z 3 1 m 2 / 

T . (s(l + /3) z G z 5 \ . ( s(l + 13) z e T z 6 \ 
L%1 ^ {—Iz^' h m 2 + LH ^ { -2^-> T ' m 2 ) 



Li m ( "A ~, 4) - ™(*. ~ 2l < - I.) ft + Mr + hh+ 

2z 3 z 5 z 6 m 2 ) \2 



I 2 I 2 I 2 I 2 

— — ItIz3 + Wz4 + — + hip + ItI/3 + lxlf3 ~ hdp + ~ + 

2U 2 (-x) + 2U 2 (x) - Li 2 + Li 2 Q + Li 2 + C(2) 

At the very end of the expression one finds an explicit imaginary part. Since the whole 
expression must be real this clearly indicates that the same imaginary contribution 
with opposite sign must be contained in multiple polylogarithms, e.g., some of them 
are sitting on branch cuts. This is in fact true for the multiple polylogarithms 

/ m 2 xz 3 —T \ ( m 2 xz 3 2T \ 

Li 3 a — — , —z- , Lisa 



sTI3 ' m 2 x) ' °' V sTI3 ' t(l - (3) J ' 

^ . / m 2 xz 3 sT[3 z 3 \ ^ . ( m 2 xz 3 sT[3 z e \ 
' ' y sT/3 ' m 2 xz 3 ' s(3 J ' ' \ s Tf3 ' m 2 xz 3 ' t(3 ) 

Indeed, one finds that the imaginary contributions cancel out when one numerically 
evaluates the result. 

As regards the length the representations of ReD^ in terms of L functions in 
|Tj and in terms of multiple polylogarithms are of similar size. The representation in 
terms of L functions contains 43 different L function expressions against 59 different 
multiple polylogarithm expressions. 
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